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We report on some recent progress made in understanding weak matrix elements of mesons in the context of the 
next-to-leading order of the large- A'^c approximation to QCD. Specifically, we first use the example of the weak 
contributions to the vr"*" — tt" mass difference to exhibit how a systematic matching can be achieved analytically 
between short distances and long distances within our large- A^c framework. We are then also able to compute 
matrix elements of the operator Qr, as they turn out to depend on the same QCD correlator as the previous pion 
mass difference. As a final example we determine the chiral counterterms governing the pseudoscalar decay into 
a lepton pair, where we briefly comment also on the special case Kl — > fJ,'^fi~ ■ This report covers the material 
presented by the authors in three separate talks at the QCD '99 conference in Montpellier. 



1. INTRODUCTION 

At the present level of accuracy, the simple de- 
scription of electroweak interactions in the Stan- 
dard Model gives an excellent fit to the experi- 
mental data at high energies collected over the 
last years, for instance at LEP ||l|. For low-energy 
phenomenology however, the description in terms 
of the fields which appear in the Standard Model 
Lagrangian is inappropriate, ft is more conve- 
nient to use an effective Lagrangian description 
where the heavy degrees of freedom of the Stan- 
dard Model are integrated out in the presence of 
the strong interactions. This procedure leads, 
however, to a rather complicated structure for 
the weak processes involving hadrons. The effec- 
tive Hamiltonian that describes the non-leptonic 
AS* = 1 decay of kaons, for instance, is given by 
a set of four-quark operators Qi, 



w"' = 5]a(Ai)g.(M), 



(1) 



modulated by the functions Ci{fi), the Wilson 
coefRcients, containing the information from the 
short-distance physics. The matching of the effec- 
tive theory to the underlying theory of the Stan- 
dard Model is performed at a high scale /Ism ^ 
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jTit ^ M\Y, where the corresponding Wilson co- 
efficients Ci{fisM) can be accurately computed in 
perturbation theory. In order to evaluate a weak 
matrix element at much lower energies fi <C (J-sm, 
the Wilson coefficients are evolved downward us- 
ing their renormalization group properties ^J . 

Physical observables evaluated with the effec- 
tive Hamiltonian in Eq. (|l|) are independent of 
the arbitrary scale ^ that was introduced in order 
to separate the short-distance physics, contained 
in the Wilson coefficients Ci{^), from the long- 
distance contributions described by the hadronic 
matrix elements < h'\Qi{^)\h > of the four- 
quark operators Qi{fJ.). However, depending on 
the scale fi at which one is working, the contri- 
bution of a given operator Qi{fi) can be more 
or less important, according to the behaviour of 
the modulating factor Ci{^) under the renormal- 
ization group evolution. The latter is available 
only within a perturbative expansion in powers 
of the strong coupling constant as ■ Although the 
present state of the art includes NNLO contribu- 
tions O] , it can only provide a reliable descrip- 
tion of the scale dependence of the Wilson co- 
efficients down to a scale fj, which may, at best, 
be taken slightly below the charm quark mass, 
/i < TTie. On the other hand, the computation of 



the matrix elements < h'\Qi{^)\h > themselves 
requires non-perturbative methods. One possi- 
bility is to have recourse to numerical compu- 
tations based on a discretized lattice version of 
QCD. Another approach consists in implementing 
in an analytic way some particular but systematic 
expansion scheme, like chiral perturbation theory 
(ChPT) and/or the large- A'^c expansion, etc. The 
fact that the final result for the matrix element 
< h'\H°"\h > must then be independent of the 
factorization scale (at the given order of the non- 
perturbative expansion considered for the eval- 
uation of the matrix elements of the four-quark 
operators) provides thus a non-trivial constraint 
and a crucial check of the whole calculation. 

Finally, at a very low scale /x <C 1 GeV, the in- 
teractions of the light pseudoscalar mesons can be 
described within a systematic expansion in pow- 
ers of momenta (chiral expansion). Indeed, in 
the chiral Hmit, these pseudoscalar states are the 
Goldstone bosons of spontaneous chiral symme- 
try breaking, so that their interaction becomes 
small at low energy. To lowest order in ChPT, 
their strong and weak matrix elements can be ob- 
tained, in the chiral limit, from the 0{p^) chiral 
Lagrangian 
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where £g accounts for the strong interactions 
between the pseudoscalars 



,(2) 
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while the weak \AS\ = 1 transitions receive con- 
tributions with the (8,1) quantum numbers of the 
chiral group SU{3)l x SU{3)b. 



,(2) 
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or with the (27,1) quantum numbers 
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Here we use the notation C^ = —i^U^D^jU , 
and as usual U denotes a unitary 3x3 matrix 



describing the octet of pseudoscalar fields. All 
strong and weak matrix elements involving only 
the Goldstone bosons can thus be evaluated at 
0{p^) accuracy and expressed in terms of only 
three low-energy constants, the decay constant of 
the pion in the chiral limit Fq, and the two con- 
stants (78 and 1727 which, because of CP violation, 
are complex numbers. For instance, in this ap- 
proximation, we obtain the following ratio of the 
K — > (tttt)/ amplitudes for the decay of kaons into 
two pions in a prescribed isospin channel / = 0, 2 
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(6) 



Higher order or quark mass corrections can, of 
course, be considered as well, but at the expense 
of introducing additional low-energy constants. 

Unfortunately, an a priori knowledge of the 
values of the constants g^ and 527, which would 
lead to a quantitative understanding of the ratio 
in Eq. (^), i.e. of the A/ = 1/2 rule, is not avail- 
able at present. Such a knowledge would require 
that one is actually able to continue the process 
of integrating out the high energy modes in (|l|) in 
a non-perturbative way down to the very low scale 
where the description in terms of the effective La- 
grangian (^ is valid. On the other hand, these 
low-energy constants can be identified in terms 
of QCD correlators of bilinear quark operators. 
For instance, Fq is exactly given by the following 
two-point correlator in the chiral limit and at zero 
momentum transfer: 



JimJ-Q2ni^(g2)) ^ F^, 
where 



(7) 



nf^(g) = 2z d'xe^'^-^{0\T{L'^ix)R''{0)^)\0) 



{q^q''-g^^q')ULR{Q' 



(8) 



with Q^ = —q^, and where L^ and R'^ are left- 
and right-handed currents with the appropriate 
flavour quantum numbers e.g. L^ = ULji_idL, 
Rp, — uiilfida. The low-energy constants of 

the strong interaction part £g' = X]n>i'^s 
describing higher orders in the chiral expansion 
can likewise be expressed in terms of the low- 
momentum behaviour of QCD correlators. For 



instance, one of the Gasser-Leutwyler constants 
J3| at order C'(p'*), Lio, is given by the coef- 
ficient of the secont term in the expansion of 
-Q^TIlr{Q^), i.e. 



Q^IiLRiQ^) = Fo + 4iioQ' + O(Q') : 



(9) 



modulo chiral logarithms that we have not dis- 
played explicitly. In general, one can find an 
analogous representation for the other low-energy 
constants of £5* as coefficients of a Taylor ex- 
pansion around zero momentum of certain Green 
functions, but not always two-point Green func- 
tions. Using dispersion relations, one may then 
express these low-energy constants in terms of ob- 
servables of the hadronic spectrum in the chiral 
limit. As a matter of fact, the low-energy con- 
stants of interest, such as Fq, and the correspond- 
ing Green's functions, such as n^'^((7), happen to 
be order parameters of chiral symmetry break- 
ing. This means that these Green's functions, 
which receive no contribution from the pertuba- 
tive QCD continuum, behave smoothly at short 
distances, and that the corresponding dispersion 
relations converge therefore rapidly. In practice, 
one may thus expect to need information on the 
hadronic spectrum only over a finite energy range, 
say up to ~1 GeV. Unfortunately such data, even 
away from the chiral limit, are not always avail- 
able. An enormous simplification comes about 
after taking the large- A^c limit of QCD. In this 
limit the spectrum consists of an infinite number 
of zero-width mesonic resonances. For instance, 
^lr{Q^) is then given by the pion pole and by an 
infinite sum over single particle vector and axial- 
vector states, i.e. 
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which, upon expanding the above large- A^c repre- 
sentation of —Q'^Hlr{Q^) around Q^ = 0, yields 
the following expression for Lio H 



4Lio = E/l-E/^ 



V ' 



(11) 
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in terms of the parameters of the zero-width 
mesonic resonances. 



The situation is more involved when it comes 
to the constants occurring in the |A5| = 1 sec- 
tor, such as 58 or 1727. Although they can still be 
expressed in terms of QCD four-point functions 
of quark bilinears which are order parameters of 
spontaneous chiral symmetry breaking, they in- 
volve for instance two left-handed charged cur- 
rents convoluted with a free W-boson propagator. 
Therefore, the whole range of momenta, and not 
only the low-momentum region, is involved. The 
same happens to the constants describing electro- 
magnetic corrections to strong processes. 

The crucial question is of course whether this 
interpolation at large Nq can be realized in a 
way that provides the correct matching of the 
scale dependence of the weak matric elements 
< h'\Qi{^)\h > with the scale dependence of the 
short-distance piece encoded in the Wilson coef- 
ficients, Ci{fj,). As we shall see, there are some 
observables where we have already been able to 
realize this program. 

2. A CLASSIC EXAMPLE REVISITED 

In the chiral limit the tt, the K and the rj 
form an octet of massless Goldstone bosons, pro- 
vided electroweak interactions are neglected. In 
the presence of electromagnetic interactions, how- 
ever, the electrically charged members of the erst- 
while octet of Goldstone bosons acquire a mass 
even in the chiral limit. This mass term can be 
described by the effective operator 
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where Qr^l are matrices of charges in flavour 
space governing the couplings of quarks to the 
photon, Qij^Qi = diag[2/3, — 1/3, — 1/3], and a 
is the fine structure constant. The constant C-y 
results from integrating out all degrees of freedom 
but the Goldstone bosons and the photon in the 
Standard Model Lagrangian. The fact that the 
QCD part is described in terms of quarks and 
gluons and not in terms of hadrons is of course 
a good part of the difficulty in performing this 
integration. 

In this classic example, it is well known that 



C~f is determined by the same two-point func- 
tion ^lr{Q^) that we have aheady encountered 
^, but contrary to Fq or to Lio, which describe 
its low-momentum behaviour, the identification 
of C~f involves an integral of ^lr{Q'^), weighted 
with the free photon propagator, over the whole 
range of momenta. 
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We stress here that Cy is thus very akin to the 
coupling constants in the chiral Lagrangian de- 
scribing electroweak interactions of hadrons, such 
as gs or (727 for instance. 

The integral in ( [l3|) can be split-off at an ar- 
bitrary intermediate scale A, but large enough 
that perturbation theory sets in 0. The high- 
momentum region of integration, Q^ ^ A^, can 
be evaluated upon taking for TlLRiQ^) its asymp- 
totic behaviour at short distances, which in QCD 
is given by (this property ensures the convergence 
ofEq. (g) 

(14) 



lim Q^IiLBiQ^) = 0, 



and 



lim Q^nLfi(Q2) = 

-47r2(^-|-0(a2))(Vi^)2 



(15) 



Notice that in the last formula, we have already 
replaced the relevant four-quark condensate by its 
factorizcd large- iVc expression. For the large- iVp 
limit representation in Eq. (lOh, these properties 
translate into the first and second Weinberg sum 
rules g 

^/^M^-^/lMl = F„^ (16) 

V A 

J2fyM^-J2flMX = 0, (17) 
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and into the relation g] 
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respectively. The corresponding contribution to 
C-y from the high-momentum part of the integral 
reads 
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For the low- momentum region of the integral, the 
large- TVc representation of Eq. ( p^ ) gives 

C<(A) = 
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after using the two Weinberg sum rules (Iq) and 
(O). Upon adding the two pieces, one obtains a 
perfect matching to the given order in 1/A, i.e. 
C<(A) -K C>(A) ~ Ci(l/A4), due to the sum-rule 
in Eq. (p^). Clearly this cancellation can be ef- 
fected to any order in 1/A, as Eq. (^ has no A 
to begin with, but the point is that, in order to 
accomplish this A matching, care must be taken 
that the constraints imposed by the OPE are sat- 
isfied to that order in 1/A [^. This lesson is very 
important in the case of weak matrix elements. 

Let us notice here that a simple ansatz that 
realizes all the properties listed above consists 
in taking into account only one vector and one 
axial- vector resonance (satisfying even some other 
stronger constraints discussed in Ref . |g] ) that we 
shall call the Lowest Meson Dominance (LMD) 
approximation to large- iVc QCD. In this case one 
has the simpler expression H 
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and 



where M^ = ^{A^F^f, Mai = V^Mp 

Fq ~ 87MeV is an estimate of the pion decay 

constant in the chiral limit HS]. Combining 



Eqs. (0|l|) and (|l]) one obtains 



2 p2 24n/6 

m_+ = Fn — z — 7ra In 2 
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I.e. TO7r+ — TItt" 



i! 5 MeV to be compared to 

4.6 MeV. 

In order to create a situation which is more 
similar to the one encountered in the case of the 
I AS" I — 1 transitions, let us rather consider the 
masses acquired by the pseudoscalar octet under 
the influence of the weak neutral currents [ |lO[ . 
They result from an effective term similar to fl2|) 



£off 



mass-.Z'^^ 



4T:aCztr{QRUQLU^) 



(23) 



but with a coupling constant Cz which involves 
the integral of I1lr{Q^) weighted by the propa- 
gator of the massive 2'°, 



^- = ii/ 'Q 



Ml 



g'nifl(Q2 



(24) 



The change of sign as compared to Eq. ( [131) is a 
consequence of the SU{2) x C/(l) quantum num- 
bers of the quarks. Using an Euclidean momen- 
tum cutoff A, one may again split the integral into 
a low-energy piece, which gives 
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and a high-energy piece, evaluated with the lead- 
ing short-distance behaviour given by Eq. (l5h, 



C|(A) 



3 1 /7,n2 1 Ml 
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Adding Eqs. (|2^) and (26) one obtains the final 
result which is, of course, A independent. 

A more common approach is to construct an 
effective Lagrangian in which the Z'^ has been in- 
tegrated out and the Green functions are renor- 
malized using the MS scheme. The relevant term 
in the Lagrangian of the Standard Model which 
is responsible for the Z-induced contribution to 
the tt"*" — tt'^ mass difference is the neutral current 



interaction term 
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(27) 



When looking for the induced effective La- 
grangian of order 0{p'^) which contributes to 
Goldstone boson masses, it is sufficient to con- 
sider left-right operators. In the absence of 
the strong interactions, the effective four-quark 
Hamiltonian which emerges after integrating out 
the Z field is 



-1 e2 
"° " Ml cos2 e 



w 



1 



sin 0wQlr ~ 2 ilLlp^T^qL) {qRl^QRqR) 



M\ 



J^LR. 



1 1 



Ml cos2 Bw 6 



X! i^t^i {qr^^Qmr) 



where 



Qlr = (qLjf^QLqL) {qRl^QRqR) , 



(28) 



(29) 



and summation over quark colour indices within 
brackets is understood. In fact, to ©(p"), only the 
first term proportional to the four-quark operator 
Qlr can contribute. In the presence of the strong 
interactions, the evolution of Qlr from the scale 
Ml down to a scale ii^ can be calculated in the 
usual way, provided this ^^ is still large enough 
for a perturbative QCD (pQCD) evaluation to be 
valid. In the leading logarithmic approximation 
in pQCD, and to leading non-trivial order in the 
l/Nc expansion, the relevant mixing in this evo- 
lution which we need to retain is simply given by 
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where fj, is the MS scale and Dul denotes the 
four-quark density-density operator 
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with Cq and Cq' the quark charges in units of the 
electric charge. This can be seen as follows: in the 
MS renormalization scheme, the full evolution of 
the Wilson coefficients cq and cd of the operators 
Qlr and Dlr. at the one-loop level is governed 
by the equations (subleading contributions in the 
1/Nc expansion have been neglected) 
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with boundary conditions: cq{Mz) — 1 and 
cd{Mz) = 1^- The resuh in eq. @ follows 
when taking cd{Mz) — 0, which is appropriate 
when keeping the one-loop leading log only, and 
from the off-diagonal term in the (transposed) 
anomalous dimension matrix. 

We are then confronted with a typical prob- 
lem of bosonization of four-quark operators. The 
bosonization of ZJ^x is only needed to leading or- 
der in the l/Nc expansion. To that order and 
to order 0{p'^) in the chiral expansion it can be 
readily obtained from the bosonization of the fac- 
torized density currents, with the result g 



DRL = ^eqeq'{q'i^qR){qRqL) 

(M)\r(t/Q.t/tQ,) 



(33) 



We find that the overall contribution of the term 
proportional to the Diil{p^) four-quark opera- 
tor, which we denote CzIdrl j is given by the ex- 
pression 



-3 1 - M^ 

CzlonL = — 77J"'^(^V')'log^ 



877 M| 



(34) 



and it is exactly the same result as the one coming 
from the d = 6 term of the OPE in the previous 



^See e.g. the lectures in Ref. |]ll[ and references therein. 



calculation of Eq. ( P6D , except for the difference 
between /i and A. 

The problem is then reduced to the bosoniza- 
tion of the operator QLfi(A'^)- We are confronted 
here with a typical calculation of a hadronic ma- 
trix element of a four-quark operator, in this case 
the matrix element (Tr+IQLfilM^)!'"'^)- The fac- 
torized component of the operator Qlr, which is 
leading in 1/Nc, cannot contribute to the 0(p°) 
term of the low-energy effective Lagrangian. The 
contribution we want from this matrix element 
is therefore the next-to- leading one in the 1/Nc 
expansion. 

The calculation proceeds along much the same 
lines as first suggested in papers by Bardeen, 
Buras and Gerard Q sometime ago ^, except 
that we shall go beyond loops generated by 
Goldstone particle interactions alone in order to 
achieve a correct matching with the logarithmic 
scale dependence of the short-distance contribu- 
tion in eq. (|3J) . We can evaluate now the matrix 
element {^^^\QLRi^J''^)\^^~^) within the framework 
of an effective Lagrangian which is a straightfor- 
ward generalization to an arbitrary number of 
massive J^ = 1~ and J^ = 1+ mesonic states 
of the effective Lagrangian corresponding to the 
LMD approximation to QCD(oo) recently dis- 
cussed in ref. ||9|]. Furthermore we shall do so 
by using an Euclidean momentum cutoff to study 
the issue of regularization. Wc leave the details 
of the calculation to our work in Ref. |1^ . The 
final result reads 
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The contribution of the Goldstone bosons alone 
corresponds to the two terms in the first line of 

^See also Refs. |l3l-|la] and references therein for more 
recent work. 



the r.h.s. of Eq. (|35|). They display a typical 
polynomial dependence with respect to the cut- 
off fj,, which can hardly provide a reasonably good 
matching with the logarithmic scale dependence 
coming from the short-distance contributions in 
CzIdrl ■ I'^ i&ct, in an MS regularization scheme, 
as commonly chosen for the evaluation of the 
short-distance Wilson coefficients, these power di- 
vergences will automatically disappear. Simply 
adding higher resonances does not by itself solve 
the problem of matching the long and the short 
distances either; however, when the information 
coming from the two Weinberg sum rules (16 17) 
and the sum rule {Kw is taken into account, the 
result of Eq. (Ba) can indeed be recast into a form 
which reproduces Eq. pa) . Notice that in an MS 
regularization scheme, the integral in eq. (35) 
should have been understood in n = 4 — e dimen- 
sions and therefore multiplied by /u'*~", /x being 
the A/5 regularization scale; and, of course, all 
the power divergences should have been put to 
zero. The result, when combined with Eq. (34) 
finally yields the contribution of Eq. (g^) plus 
Eq. (|2^) as before. We insist on the fact that, 
regardless of the regularization one chooses, the 
calculation we have done of the matrix element 
(7r^|QLi?(^^)|7r+) is an exact calculation to next- 
to-leading order in the l/Nc expansion and to 
0(jP) in the chiral expansion. 

In Fig. 1 we compare the complete result (with 
just one resonance for the vector and axial-vector 
channels) with the value obtained for Cz when 
only Goldstones arc included in the long-distance 
part, i.e. only the first two terms of Eq. (35) 
are kept. One immediately sees that the com- 
plete result shows a flat dependence on the cutoff 
/J,, as it should, while the result with only Gold- 
stones shows a certain dependence on fi. One 
also sees that for a reasonable value of the cutoff 
H ^ 1.15 GeV (which is some sort of average be- 
tween the p and oi masses) the result with only 
Goldstones gives the right answer but that, if one 
misses the right value of /i in this guess, one can 
easily even flip the sign of the whole contribution. 

We think this clearly exemplifies the potential 
dangers of a calculation with only Goldstones and 
a physical Euclidean cutoff |13] [15| . 
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Figure 1. Plot of the total contribution to Cz, Eqs. 
d34) plus ( 351) , i) when only Goldstones are considered 
in Eq. (BS) (dashed curve), ii) full result including all 
terms in Eq. ( |3q ) with just one vector and axial- 
vector resonances (solid curve). 



3. ELECTROWEAK PENGUIN OPERA- 
TORS 

Within the framework discussed above, we have 
also shown n& that the K — > tttt matrix elements 
of the four-quark operator. 
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can be calculated to first non-trivial order in the 
chiral expansion and in the 1/Nc expansion. 

The operator Q^ emerges at the Mw scale from 
considering the so-called electroweak penguin dia- 
grams. In the presence of the strong interactions, 
the renormalization group evolution of Qj from 
the scale Mw down to a scale /i < rric mixes this 
operator with others, in particular with the four- 
quark density-density operator Qg 



= -12 



X! eq{sLqR){qRdi 



(37) 



These two operators, times their corresponding 
Wilson coefficients, contribute to the lowest order 
0{p'~') effective chiral Lagrangian which induces 
jAS*! — 1 transitions in the presence of electro- 
magnetic interactions to order 0{a) and of vir- 



tual Z'^ exchange, i.e., the Lagrangian 111 
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where A}^ is the effective left-handed flavour 
matrix ( A 



(23) \ ^ 



St2S3j 



is the only possible invariant which in the Stan- 
dard Model can generate \AS\ = 1 transitions to 
orders 0{a) and 0{p'^) in the chiral expansion. 
The coupling constant h is dimensionless and, a 
priori, of order 0{N^) in the I/Nq expansion. 
It plays a crucial role in the phenonienological 
analysis of radiative corrections to the K — > tttt 
amplitudes, hence the interest of identifying all 
the possible contributions to this constant. 

The bosonization of the operator Qj to next- 
to- leading order in the 1/Nc expansion turns out 
to be entirely analogous to the bosonization of 
the operator Qlr = {qLYQhqL) {qRj^VRqn.) 
which governs the electroweak 7r+ — tt*^ mass dif- 
ference discussed in the previous section. Because 
of the LR structure, the factorized component of 
Qt, which is leading in 1/Nc, cannot contribute 
to the 0{p^) low-energy effective Lagrangian in 
Eq. (p8|). The first 0{p^) contribution from this 
operator is next-to-leading in the 1/Nc expansion 
and is given by the integral |16| , 

d'^q 



Q 



'Sigfiu 



(27r)4 

(23) 



Wr(.i) X 



tr(ux'^Pu^QR 



(39) 



involving the same two-point function as in 
Eq. (ra). This integral, however, is divergent for 
large Q^ and needs to be regulated. The usual 
prescription [0| for the evaluation of integrals 
such as this, consists in taking a sharp cut-off 
in the (Euclidean) integration over Q^, 



Qi 



-6- 
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dg^O^ (_g2n^^(g2)) 



tr[/Af^ C/tQfl .(40) 



Inserting the same large-iVc expression for the 
function ^lr{Q^) as in Eq. (p^, with the short- 
distance constraints between the couplings and 



masses of the narrow states incorporated, the in- 
tegral on the r.h.s. of Eq. (|4(]| ) becomes then only 
logarithmically dependent on the ultraviolet scale 
A, and one obtains the following result 



(z,j = 1,2,3). This ^/iMllog 



dQ'Q- 

Ml 



-Q'^LRkO')) 



E/'*^'i°g 



A2 
Ml 



(41) 



for values of the cutoff much larger than any res- 
onance mass included in the difference. Notice 
that if only the contribution from the Goldstone 
pole had been taken into account, the resulting 
expression would have displayed a polynomial de- 
pendence on the cut-off scale A. Another pos- 
sibility is to evaluate the integral in Eq. (M) 
within a dimensional regularization scheme, say 
MS, in which case one obtains the same result 
as in Eq. (|4l|), but with the correspondence be- 
tween the cut-off A and the MS subtraction scale 
/i given by A = ^ • ee . In this case we can check 
that our result satisfies the correct renormaliza- 
tion group equation, i.e. by acting with [idjdyi, on 
Q'j in Eqs. (|0 41) one obtains 






dfi 




(42) 



where one has to bosonize the operator Qg in Eq. 
(37) utilizing, mutatis mutandis, the rule of Eq. 
(33) and the constraint of Eq. (|l8|). 

The bosonic expression of Qy given by Eqs. ([40| ) 
and (^) enables us to compute the _ftr — > tttt ma- 
trix elements induced by this operator which, fol- 
lowing the usual conventions, we express in terms 
of the following isospin amplitudes 

((7r7r),|Q7|i^°), 



(Q 



7)1 



0,2. 



(43) 



To leading order 0(jP) in the chiral expansion 
and to next-to-leading order in the 1/A^c expan- 
sion, 0{l/^/Nc) for K ^ TTTT amplitudes, we ob- 
tain the result 

(g7>o = V2(g.), = ^M^x 




Y^fyM' 




,(44) 



It has become customary to parameterize the 
resuhs of weak matrix elements of four-quark op- 
erators Qi in terms of the factorized contribu- 
tions from the so-called vacuum saturation ap- 
proximation, modulated by correction factors, the 
so called B-factors. Although the resulting B fac- 
tors for the A/ = 1/2 and A/ — 3/2 transitions 
generated by the Q7 operator are found to de- 
pend only logarithmically on the matching scale 
fi, their actual numerical values turn out to be 
rather sensitive to the precise choice of /i in the 
GeV region. Furthermore, because of the normal- 
ization to the vacuum saturation approximation 
inherent to the (rather disgraceful) conventional 
definition of B-factors, there appears a spurious 
dependence on the light quark masses as well. In 
Fig. 2 we show our prediction for the ratio 



n(3/2) ^ \" '-^' I"' / (AK') 

versus the matching scale /i defined in the MS 
scheme. This is also the ratio considered in recent 
lattice QCD calculations |lq|. [In fact, the lattice 
definition of B-^ uses a current algebra relation 
between the K — > tttt and the K —^ n matrix 
elements which is only valid at order 0(p°) in 
the chiral expansion.] In Eq. (Eq), the matrix 
element in the denominator is evaluated in the 
chiral limit, as indicated by the subscript . 

It would be much better, whenever possible, to 
compare lattice results directly of the amplitude 



\Q\ \K^) with our prediction 



(^+ig 



(3/2) 



K+) = 






yflM%\og 



Ml 



A2 



M2 



A2 



M2 



Ml 



(46) 



where in the second line we have used the LMD 
approximation discussed in Ref . Q , or the equiva- 
lent expression in the MS regularization, namely. 
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Figure 2. The B) factor in Eq. (Es) versus ^ in 
GeV. Solid lines correspond to (ma-|-7nd)(2 GeV) = 
158 MeV; dashed lines to (ms -I- md)(2GeV) = 
100 MeV. 



3A4 
47r2 



2 log 2 -log 



Ml 



(47) 



Numerically one obtains —0.02 ± 0.01 GeV^ for 
the expression in Eq. (^^ evaluated at /i = 2 
GeV. The error is an estimate of 1/A^c corrections 
and corrections to the LMD limit. 

It is interesting to compare the A dependence 
of Eq. (H^) with the ^ dependence of Eq. (|7|). 
The two results clearly coincide for asymptotic 
values of their respective A and /z scales. The 
situation, however, is rather different for values 
of these scales in the GeV region. This can be 
best seen by looking at the functional relation 
between these two scales which follows from iden- 
tifying the two expressions in Eqs. ( pq ) and (p^: 
setting X = ■^rre^''^ and v = ttt this relation is 

4 illy/Hi" 



given by the function x 



i+y 



The require- 



ment that y >0 results in a non trivial constraint 
X > 4, i.e., /i > 2e-^/^Mv = 1.69My. In fact, 
at the value /i = 1.69My the matrix element in 
Eq. ( [47I ) flips its sign in contradiction with a gen- 
eral positivity property |19|| which demands that 
-Q2n(Q2) > for aU values of < g^ < 00. Al- 
though the specific critical value /i — 1.69Afv de- 
pends on the hadronic LMD approximation which 
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we have made, it shows that pushing the match- 
ing of four-quark operators at low /x^ values in 
the MS regularization scheme may be very dan- 
gerous and can lead to totally misleading results. 

4. DECAY OF PSEUDOSCALARS INTO 
LEPTON PAIRS 

The examples discussed so far involved only 
the two-point function 11/,^. In this section we 
present an example which allows us to study the 
matching between short and long distances in 
the case where a three-point function is involved 
]20[ . The physical processes of interest are the 
decay of 7r° or -q into a pair of charged leptons. 
These processes are dominated by the exchange 
of two virtual photons, as shown in Fig. 3, and 
it is therefore phenomenologically useful to con- 
sider the branching ratios normalized to the two- 
photon branching ratio [P = tt'^^t]) 



R{P ^ tt 



I ttMp 



_ Br{P ^£+£-) 
' ~ Br{P -^ 77) 

J 



(48) 



with Pi [s) = y/l — 4m'j/s. The unknown dynam- 
ics is then contained in the amplitude A{Mp). To 
lowest order in the chiral expansion the contribu- 
tion to this amplitude arises from the two graphs 
of Fig. 3 with the result 



A{.s) = xp(m) + ^ 



Nc 
3 



3 M 

2 U' 



C{s) 



where 



XttO = Xv 



(Xi +X2) _ 



X 



(49) 



(50) 



with xi and X2 the couplings of the two coun- 
terterms which describe the direct interactions of 
pseudoscalar mesons with lepton pairs to lowest 
order in the chiral expansion pi[ 

X [xMQrQrD^UU^ - QlQlD^U^U) 

+ X2tr (U^QrD^UQl - UQlD^U^Qr)] . (51) 



Ct 




Figure 3. The lowest order contributions to the P — > 
l^£~ decay amplitude. The second graph denotes 
the contribution from the counterterm Lagrangian of 



Eq. {^. 



The function C{s) in Eq. ( [49| ) corresponds to 
a finite three-point loop integral which can be 
expressed in terms of the dilogarithm function 
Li2(x) = - J^{dt/t)\ii{l -t). For s < 0, its 
expression reads 



C{s) = 



1 



lids) 



Li 



4 \Hs 



Pi{s) - 1 

Pt{s) + l 

Pi{s) - 1 



1 



(52) 



The corresponding expression for s > Amf is ob- 
tained by analytic continuation, using the usual 
ie prescription. The loop diagram of Fig. 3 orig- 
inates from the usual coupling of the light pseu- 
doscalar mesons to a photon pair given by the 
well-known Wess-Zumino anomaly |22] . The di- 
vergence associated with this diagram has been 
renormalized within the MS minimal subtraction 
scheme of dimensional regularization. The loga- 
rithmic dependence on the renormalization scale 
/i displayed in the above expression is compen- 
sated by the scale dependence of the combina- 
tion x(/x) of renormalized low-energy constants 
defined above. Let us stress here that, as shown 
explicitly in Eq. (^9|) and in contrast with the 
usual situation in the purely mesonic sector, this 
scale dependence is not suppressed in the large- 
Nc limit, since it does not arise from meson loops. 
As a first step towards its subsequent evalua- 
tion we shall identify the coupling constant x in 
terms of a QCD correlation function. For that 
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purpose, consider the matrix element of the light 
quark isovector pseudoscalar density P^{x) — 
^{uij5U—dij5d)(x) between leptonic states in the 
chiral limit. In the absence of weak interactions, 
and to lowest non-trivial order in the fine struc- 
ture constant, this matrix element is given by the 
integral 

<r{p')\p'{o)\£-{p)> 

(2^)4 [{p'-q)^-mj]q^p'-p-q)^ 



X I 



-74„, Jq-x^i{p'-p-q)-y 



d'ye'' 



X <0\T{f;^{x)fr{y)P'm\0>, (53) 

with j"" = g (2m7^m — d'f/j^d — sjf^s). In the chi- 
ral limit, the QCD three-point correlator appear- 
ing in this expression is again an order parameter 
of spontaneous chiral symmetry breaking. Bose 
symmetry and parity conservation of the strong 
interactions yield 



d^a 



d ye 



x<Q\T{j;^\x)jr{y)PHO)}\0> 
2 



■eA"'a/3'?"92 ^(9?,92'(9l +92) 



(54) 



with n{qlql (gi + 92)') = ^(92', 9?, (91 + q2?)- 
For very large (Euclidian) momenta, the leading 
short-distance behaviour of this correlation func- 
tion is given by 

lim H{{Xqi)MXq2)M\qi + \q2f) 



1 

2A4 



< tjjijj > 



9? + 91 + (91 +92)' 



qlqliqi 



q2f 



' A4'A6 



(55) 



Actually, what matters for the convergence of the 
integral in Eq. ( ^3|) is the leading short-distance 
singularity of the T— product of the two electro- 
magnetic currents, which corresponds to 

lim n{{Xq)\{p'-p-\q)Mp'~pf) 
A-^ q^{p' — pY 



O 



A2' A3 



(56) 



Thus, the loop integral in Eq. (^ is indeed con- 
vergent. The QCD corrections of order 0{as) in 
Eqs. (g^ and (|5^) will not be considered here. 
Let us however notice that since the pseudoscalar 
density P^{x) and the single-flavour < ^^ > con- 
densate share the same anomalous dimension, the 
power-like fall-off displayed by Eqs. ( pq ) and ( [Sq ) 
is canonical, i.e. it is not modified by powers of 
logarithms of the momenta. 

At very low momentum transfers, the same cor- 
relator can be computed within Chiral Perturba- 
tion Theory (ChPT). At lowest order, it is satu- 
rated by the pion-pole contribution, given by the 
anomalous coupling of a neutral pion, emitted by 
the pseudoscalar source P'^(O), to the two electro- 
magnetic currents, i.e. 



H(0, 0,(91+92)') 



Nc < ijjip > 
8^ F§ 



1 



(91 



92)2 

(57) 



where the ellipsis stands for higher orders in the 
low-momentum expansion. The matrix element 
< i~{p') I p3(o) I l~{p) > itself may also be eval- 
uated in ChPT. At lowest order, it is given by 
the diagrams of Fig. 3, where the (off-shell) pion 
is now emitted by the pseudoscalar source P^{0). 
The result reads, with t = {p' — p)^, 



<r{p')\pHo)\r{p)> 



32TrH 



f2 
^0 



ChPT 

■mtu{p')-f^u{p) A{t) 



(58) 



with the function A{t) defined in Eqs. (H3) and 
([5^). The contribution of the loop diagram of 
Fig. 3 is obtained upon replacing, in Eq. (bSh, 
the three-point QCD correlator by its lowest or- 
der chiral expression given in Eq. (|5^). The cou- 
pling constant x(/i) is thus given by the residue 
of the pole at i = of the matrix clement 
< i-{p') I P3(0) I i~{p) >, after subtraction of the 
contribution of the two-photon loop, i.e. 

X(m) < -0^ > -f ,. . ^ 

— Zi m — rniu[p )"/5u{p) 



327r4 



lim 



2i 



uip'hxJrMp) X 



d'^q 



{p'-pf 



p)2^o J {2ttY [{p' - qY - ■mj]q'^ [p' - P - qY 



12 



{p' -p- 2q)a 



q"{p'-p)^- ip'-pTq' 



n{qMp'-p-qf,{p'-pf) 



H{Q,{),{p'~pf) 



(59) 



In the large- iVc limit, the three-point correlator 
( |54| ) is given by the tree-level exchanges of vector 
and pseudoscalar resonances, as shown in Fig. 4, 
so that the singularities in each channel consist of 
a succession of simple poles. This involves cou- 
plings of the resonances among themselves and to 
the external sources which, just like the masses 
of the resonances themselves, cannot be fixed in 
the absence of an explicit solution of QCD in the 
laige-Nc limit. Nevertheless, the general struc- 
ture of the quantity appearing in Eq. ( p9| ) is of 
the form 




7 em 7 em 




V.P 






) 



Figure 4. The contributions to the vector-vector- 
pseudoscalar three-point function in the large-A^c 
limit of QCD. The sum extends over the infinite num- 
ber of zero-width vector (V) and pseudoscalar (P) 
states. 



hm ip'~pfniq^q'Ap'-pf)^ 



Y^M, 



1 < tpip > 

ay 
[q^ - M^) ~ [q^ M^Y 



bvq' 



, (60) 



where a priori the sum extends over the infi- 
nite spectrum of vector resonances of QCD in 
the large- A^c limit. Equation (pOh follows from 
the fact that its left-hand side enjoys some addi- 
tional properties: i) In the pseudoscalar channel, 
only the pion pole survives, while massive pseu- 
doscalar resonances cannot contribute, ii) The 
momentum transfer in the two vector channels is 
the same, iii) Its high-energy behaviour is fixed 
by Eq. (H). 

Even though the constants ay and by depend on 
the masses and couplings of the vector resonances 
in an unknown manner, they are however con- 
strained by the two conditions 



E 

V 



ay 



Nc 
47r2' 



^(av 



^Mt 



2Fl (61) 



which follow from Eqs. (57) and (|56|), respec- 
tively. Notice that there are no contributions 
from the pcrturbative QCD continuum to these 



sums. Taking the first of these conditions (which, 
coming from the anomaly, has no 0{as) correc- 
tions) into account, we obtain 



x(m) 



57V£ 
12 



-27r 



^E 

V 



1 \^^v 
ay In — — 



-bx 



.(62) 



This equation, together with the two conditions 
(|6l|), constitutes the central result of Ref. p^ . 
This is as far as the large- A^c" limit allows us to 
go. Let us point out that the scale dependence 
of x{^J') is correctly reproduced by the expression 
( |6^ ) , again as a consequence of the first condition 
in Eq. (IH). 

Within the LMD approximation of large-iVc 
QCD, it is easy to write down the expression 
of the correlation function H((/f,g|,(gi -I- 92)^) 
which correctly interpolates between the high en- 
ergy behaviour in Eq. (|2^) and the ChPT result 
in Eq. (1^) II 



/LMD/ 



92 j 



W""^(gf, 92^(91 

(gi + 92)^ 



ql 



ql 



= -2 < V'V' > 

M4aLMD/^^2 



{ql-Ml){ql 



Ml){qi 



92 j 



(63) 



Notice that this expression also correctly repro- 
duces the behaviour in Eq. ( p6[) . In this approx- 
imation, the two conditions (|6l|) completely pin 
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down the two quantities ay and by in terms of 
Fq and of the mass My of this lowest lying vec- 
tor meson octet, 



,LMD 



N, 



N, 



2F-^ 



and by - 2 • 104J 



47r2 



47r2 Af2 



With the results of Eq. {<§§, and for 7Vc=3, it 
follows from Eq. (H) that 



11 3^^ (Ml 
A* 



4,'^. (65) 



Numerically, using the physical values Fq : 
87 MeV and My = A-h = 770 MeV, we obtain 



X^^^^^{p = My) = 2.2 ±0.9, 



(66) 



where we have allowed for a systematic theoreti- 
cal error of 40%, as a rule of thumb estimate of the 
uncertainties attached to the large- iVc and LMD 
approximations. The predicted ratios of branch- 
ing ratios in Eq. (^) which follow from this result 
|26[ are displayed in Table 1. We conclude that, 
within errors, the LMD-approximation to large- 
Nc QCD reproduces well the observed rates of 
pseudoscalar mesons decaying into lepton pairs. 

At first sight this approximation may resemble 
good old Vector Meson Dominance (VMD). How- 
ever, there are important differences. Firstly, the 
systematic use of the 1/Nc expansion in QCD 
justifies the use of single particle intermediate 
states (and substantiates what otherwise is only 
an ansatz in VMD). And secondly, the use of the 
OPE resolves certain difficulties in the traditional 
VMD phenomenological approach, such as e.g. 
the ambiguity in the use of the VMD prescrip- 
tion for just one photon or for both photons in 
the decay P —f t^(.~ . This ambiguity is impor- 
tant, since making the wrong choice may result in 



Table 1 

The values for the ratios R{P -^ l^ (.~) obtained 
within the LMD approximation to large-A^c QCD 
and the comparison with available experimental re- 
sults. 



R 



LMD 



~R{^ 

R{v- 



Experiment 



->e+e~) X 10" 

M+At") X lO'^ 
>e+e-) X 10** 



6.2 ±0.3 7.13 ±0.55 |27| 
1.4 ±0.2 1.48 ±0.22 |6| 
1.15 ±0.05 — 



a violation of the OPE constraints given by Eqs. 

(HSii- 

The situation in the case of the \AS\ = 1 de- 
cay K1 -^ £^i~ is slightly more delicate. It has 
recently been shown [^ that these processes can 
also be described by the expressions (^) and (^9|) , 
but with an effective constant xj^o containing an 
additional piece from the short-distance contri- 
butions [0. Of course, a cast-iron understand- 
ing of these transitions is very important |32| 
as the evaluation of xM could then have a po- 
tential impact on possible constraints on physics 
beyond the Standard Model. At present, the 
most accurate experimental determination of the 
K^ -^ ^Ji^ pT branching ratio p3] gives the result: 
Br{Kl -> /i+/x-) = (7.18 ±0.17) x lO'^. Using 
the experimental branching ratio pq| Br{K^ — > 
jj) = (5.92 ±0.15) X lO"'', this leads to a unique 
solution for an effective Xk° — 5.17 ± 1.13. The 
authors of Ref. H^] argue that, to a good approx- 
imation, one has 



Xk° = X - 



^77 



SxsD , (67) 

where x is the constant defined in Eq. (^0|) and 
e-y^ is the sign of the on-shell form factor c(0,0) 
(for the notation, we refer to BO]; see in particu- 
lar Eqs. (6) and (7) of that reference) describing 
the K^ -^ 77 decay. However, as also discussed 
in Ref. Ig^l i the decay Kl -^ 77 itself is rather 
problematic. The problem comes from the fact 
that at lowest order in the SU{3)l x SU{3)r chi- 
ral expansion, one obtains c(0, 0) — 0. There are 
two ways to bypass this situation, either by con- 
sidering higher orders, or by including the r]' as 
an explicit degree of freedom already at leading 
order, which can be done within the framework 
of the U{3)l X U{3)r combined chiral and large- 
Nc expansions along the lines suggested in Refs. 
Pq ]. In both cases, additional contributions to 
the counterterm Lagrangian in (|5lj), and thus to 
the effective constant Xk° : have to be consid- 
ered: quark mass corrections in the first case, 
1/Nc corrections in the second. To the best of 
our knowledge, an accurate estimate of either of 
these corrections has, unfortunately, not been at- 
tempted so far. The analysis in the U{3)l x C/(3)_r 
framework performed in Ref. pQ] obtains the 
value of c(0, 0) from the experimental number for 
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Kl — > 77 and determines its sign, on combined 
large-iVc and phenomenological grounds, to be 
positive. In a large-iVc calculation this experi- 
mental value for c(0, 0) would require going be- 
yond the leading term because in the strict limit 
Nc ^ oo — keeping vfiu — rud — 0, rus ^ for 
simplicity — one finds again c(0,0) = 0. Con- 
sistency demands, then, that subleading terms 
be also included in the counterterms which con- 
tribute to Xk" : something that goes beyond the 
scope of this work |g^. If, on the other hand, 
one accepts the plausible phenomenological ar- 
guments of Ref. |3y] whereby these subleading 
terms are neglected, our result (^) leads then 
to Xk'> = 0.4 ± 1.1, corresponding to a ratio 
RiKl'^-^ A^+M") = (2.24 ±0.41) x 10"^ which 
is 2.5a above the experimental value R{K^ -^ 
fi+fi-) = (1.21 ± 0.04) X 10-^ 
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